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ABSTRACT 

The average kinetic energy of the heavy quark inside B or D meson is computed by means 
of the instantaneous Bethe-Salpeter method. We first solve the relativistic Salpeter equation 
and obtain the relativistic wave function and mass of 0~ state, then we use the relativistic 
wave function to calculate the average kinetic energy of the heavy quark inside heavy meson 
of 0~ state. We find that the relativistic corrections to the average kinetic energy of the 
heavy quark inside B or D meson are quite large and cannot be ignored. We estimate 
^1{=-X^) ^ 0.35(B°,5±), 0.28 {D°,D^), 0.43(5,), 0.34(1),), 0.96 (5^) and 0.62 (t?^) 
GeV2. 
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1 Introduction 



In recent years, the study of hadronic processes involving heavy quarks has attracted continuous interest 
both in experiment and in theory. The difficulty of full theory of QCD, which is dynamic theory describing 
the quark and gluon, lead us to the theoretical achievements of the Heavy Quark Effective Theory (HQET) 
The latter describes the dynamics of heavy hadrons, i.e. hadrons containing a heavy quark Q, when 
nig — > oo. The theory is based upon an effective lagrangian written in terms of effective fields, which is a 
systematic expansion in the inverse powers of the heavy quark mass ttIq. The C'^^;^^ lagrangian reads 
as follows: 

1 - n - rr C^'^ / 1 \ 

C = Kiv DK + - — K[{iD±f]K + i^K^^^K + 0[ — ) , (1) 

where the velocity-dependent field is the heavy quark field, and is the heavy quark four-velocity 
within the heavy hadron. Then the total momentum is written as Pg = rrigV + where the residual 
momentum q is the difference between the total momentum and the mechanical momentum; = 
9^ — igA^^ is the covariant derivative, and D^'^ = — v^v ■ D contains its components perpendicular to 
the hadron velocity. In the hadron's rest frame we have (i-Dj.)^ = . The second operator appearing in 
Eq. corresponds to the kinetic energy resulting from the residual motion of the heavy quark, and the 
third one in Eq. ^ the Pauli chromomagnetic interaction operator which describes the interaction of 
the heavy quark spin with the chromomagnetic gluon field. Their matrix elements can be parameterized 
as follows P]: 

f^AHQ) = ^ , (2) 

t^oiHg) = ^ , (3) 

where Hq denotes generically a hadron containing the heavy quark Q with the usual normalization 
< Hq\KK\Hq >= 2Mh. 

These two quantities are interesting for several reasons. In the HQET, heavy hadron mass is expected 
to scale with the heavy quark mass rrig as: 

M^ = m,+A+4^ + -, (4) 
where A represents the difference between the mass of the hadron and that of the heavy quark in the 
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Table 1: Theoretical estimates of the parameter /^^ of B^^d (QCDSR: QCD sum rules, HQSR: heavy-quark 
sum rules, Exp.: experimental data on inclusive decays, QM: quark models.) 



Reference 


Method 


A [GeV2] 


H , 1 Pi'i' gI^"\7" ^Mll TTTQ l/Ll 

H/lcLoivy, OliUljclK 1^1 




u. xo zc u.uu 


Ball Braun 1^ 


OCDSR 




Neubert 


QCDSR 


0.10 zb 0.05 


Gimenez et cd. [7] 


Lattice 


n no 4-0 1/1 
— u.uy ± u.i4 


A S Kronfeld et al. IHj 


Lattice 


0.45 ± 0.12 


Bigi et al. [2] 


HQSR 


> 0.36 


Gremm et al. 


Exp. 


0.19 ±0.10 


Falk etal. ^ 


Exp. 


0.1 ^ 0.16 


Chernyak JJ| 


Exp. 


0.14 ±0.03 


M Battagha et al. [E] 


Exp. 


0.17 


Hwang et al. 


QM 


0.4 ^ 0.6 


De Fazio [H] 


QM 


0.66 ±0.13 


S Simula US] 


QM 


-0.089 


T Matsuki et al. TO 


QM 


0.238 



uiq — > oo limit. In this limit, it can be related to the trace anomaly of QCD [3]: 



2Mh 4a^ 

where (3 is the Gell-Mann-Low function. Moreover, if the inclusive semileptonic width of a heavy hadron 
is calculated by an expansion in the powers of ^— , the following results are found: the leading term of the 
expansion coincides with the free quark decay rate (spectator model); no corrections of order appear 
in the rate, and the corrections depend on /i^ and /i^ jTT]. As a consequence, these parameters enter 
in the ratio of hadron lifetimes and in the lepton spectrum in inclusive transitions, which in principle are 
quantities directly comparable with experimental data. Many authors have given theoretical estimates 
of ^7r^ and /iG^ using different methods, but different results are obtained for the estimation of (see 
Table 1). Even though there may be different definitions of these two quantities, our knowledge of them 
is still far from clear due to large discrepancies, and a more careful study is still needed. 

In this letter, we give a relativistically calculated version of ^J^, i.e. we calculate the average kinetic 
energy of the heavy quark inside heavy meson in 0^ state by means of the Bethe-Salpeter method |18j . 
We solve the relativistic Salpeter equation [19] in Section 2, and give the mass and relativistic wave 
functions of heavy meson in 0~ state in Section 3. Finally, we use these relativistic wave functions to 



2 



calculate the average kinetic energy of the heavy quark in Section 4. Discussions and conclusions are also 
in Section 4. 

2 Instantaneous Bethe-Salpeter Method 

It has been known that the Bethe-Salpeter (BS) equation is one of the frameworks to describe bound 
state systems relativistically and has a very solid basis in quantum field theory. So it is very often used 
to describe bound state problems, and even in the current literature many authors would like to base the 
constituent quark model on the BS equation. For instance, in the constituent quark model the mesons, 
corresponding quark-antiquark bound states, are described by the BS equation as: 

j^:^V{p,k,q)x{k) , (5) 

where x(g) is the BS wave function with the total momentum p and relative momentum q, and V{p, k, q) 
is the kernel between the quarks in the bound state. The momenta PqiPq are those of the constituent 
quarks 1 and 2: For a heavy meson with a heavy and a light valence quark, we can treat one of these 
two constituents as a heavy quark and the other as a light quark, e.g. we treat the quark as the heavy 
quark pi = p^ and the anti-quark as the light quark p2 = Pq. The total momentum p and the relative 
momentum q are defined as: 

Po = «iP + q, ai = 



- - - _|_ 



Pq = 012P - q, OL2 



rriQ + nig 

One can see that these definitions are just the same as in the HQET, where aip is the mechanical 
momentum of the heavy quark which describes the heavy quark moving together with the meson, and the 
relative momentum q is nothing but the residual momentum of the heavy quark inside meson. However, 
our method is not the HQET and we do not have the limit of oo, so the light quark momentum 

have the meaning analogous to that of the heavy quark. 

The BS wave function x{q) should satisfy the following normalization condition: 



-Tr 



X(fe)^ [S^HpQ)S^\pq)SHk-q) + Vip,k,q)] x{q) 



= 2^Po , (6) 



(2vr)4 

where Si{pq) and £'2(^5) are the propagators of the two constituents. In many applications, the kernel 
of the four-dimensional BS equation is "instantaneous", i.e. in the center of mass frame of the concerned 
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bound state {P= 0), the kernel V{p, k,q) of the BS equation takes the simple form: 

V{j>,k,q) V{k,q) = F(| I, I Q I, cose) , 

where is the angle between the vectors k and 1 . Then the BS equation may be reduced to a three- 
dimensional one. Compared with the conditions to solve a three-dimensional equation, i.e. to evaluate its 
eigenvalues and eigenfunctions, the conditions to solve a four-dimensional one are much more complicated. 
Thus if the kernel of the BS equation for the considered problem is instantaneous, then we always would 
like to do the 'reduction' from four-dimensional to three-dimensional. Salpeter was the first to do this 
reduction, so the reduced BS equation with instantaneous kernel is also called the Salpeter equation. Here 
we briefly repeat his method and solve the full Salpeter equation. This equation is relativistic although 
it has an instantaneous kernel, so we will obtain the relativistic wave function of bound state. 

Since in the HQET the heavy quark momentum is described by using the covariant derivative = 
dfj, — igA^, and the kinetic energy of the residual motion of the heavy quark by using a covariant form 
Dx, it is convenient to write the BS equation in a covariant form. To do this, we divide the relative 
momentum q into two parts, q\\ and q_\_, a parallel part and an orthogonal one to the total momentum of 
the bound state, respectively, 

ql^(p.q/Ml)p^^ , 
Correspondingly, we have two Lorentz invariant variables: 

In the center of mass frame P= 0, they turn out to be the usual component qo and \ q \, respectively. One 
can see that in the rest frame of bound state the orthogonal residual momentum of the heavy quark is 
just the orthogonal relative momentum, i.e. iD = q. Now the volume element of the relative momentum 
k can be written in an invariant form: 

d^fc = dkpk'^dk^dsdcp , (7) 

where (j) is the azimuthal angle, s = {kpqp — k ■ q)/{kj,qj.). The instantaneous interaction kernel can be 
rewritten as: 

V{\k-'q\)=V{k^,q^,s). (8) 
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Let us introduce the notations (fpiq'^) and T]{q'^) for three dimensional wave function as follows: 



viQ±) 



(9) 



Then the BS equation can be rewritten as: 

xiq\\,q±) = SiiPQ)viq±)S2iPg) ■ 

The propagators of the two constituents can be decomposed as: 



(10) 



Si{pi 



+ 



with 



J{i)qp + OiiMn - tOip + ie J{i)qp + aiJVIjj + uj^p - ie ' 



2u)ip 



Mh 



uJip ± J{i){mi + 



(11) 



:i2) 



where i = 1,2 for heavy quark and light anti-quark, respectively, = ujq, uj2p = Wq, and J{i) = (— l)*"*"^. 



Here Ajp(g_|_) satisfy the relations: 



tp \ 



Due to these equations, may be considered as j)— projection operators, while in the rest frame = 
they turn to be the energy projection operator. 



Introducing the notations '^^^{qA_) as: 



^^H^^) ^ K%{q^)^^p{q^)^K%{q^) , 
and taking into account j^j^ = 1, we have 

v^piq±) = ^p^iq±) + Vp'{q±) + fp^iq±) + Vp'iq±) 

With contour integration over qp on both sides of Eq. IjlOp , we obtain: 

, _ A|p(g±)??p(g±)A^p(g±) Kpiq±)r]piq±)A2piQ^) 

and we may decompose it further into four equations as follows: 

{Mh -^q- ^q)^p^{q±) = At(g±)r/p(gx)At(g±) , 



(14) 



^;-{qi_) = ^-+{qi_)=Q. (15) 

In Ref. IT^, Salpeter considered the factor {Mh — uJq — ujq) being small, so he kept the first of Eqs. 
H15|) only. It is the 'original' instantaneous approximation proposed by Salpeter and followed by many 
authors in the literature. Whereas in this paper we re-examine the BS equation with an instantaneous 
kernel, i.e. we try to deal with it exactly including the second of Eqs. (|15() . The complete normalization 
condition (keeping all the four components appearing in Eqs. for BS equation turns out to be: 



(2vr)2 



2p, . (16) 



Mh Mh Mh Mh _ 

To solve the eigenvalue equation, one has to choose a definite kernel of the quark and anti-quark in the 
bound state. As usual we choose the Cornell potential, a linear scalar interaction (confinement one) plus 
a vector interaction (single gluon exchange): 

I{r) = Vs{r) + + 7o ® 7°K(r) = Ar + - 7o ® 7°^— ' (1^) 

where A is the string constant, as{r) is the running coupling constant. Usually, in order to fit the data 
of heavy quarkonia, a constant Vq is often added to the scalar confining potential. 

It is clear that there exists infrared divergence in the Coulomb-like potential. In order to avoid it, we 
introduce a factor e~"^: 

14(r) = -(l-e-n , 
a 

A fy 

K(r) = ---e-"' . (18) 
6 r 

It is easy to show that when ar <C 1, the potential becomes identical with the original one. In the 
momentum space and the rest frame of the bound state, the potential reads: 

1(9) = 14(9) +7o®7°K(^) , 

Vs{q) = -{- + Vo)S^{q)+ ^ ^ 



K(^') = -^^f^. (19) 



(Q +a2) 



The coupling constant as{Q) is running: 



127r 



1 



log (a + ^ 



QCD 



Here the constants A, a, a, Vq and Aqcd are the parameters that characterize the potential. 

3 Heavy Mesons in State 



Following the method jSnj, the general form for the relativistic Salpeter wave function of the bound state 
= 0^ can be written as (in the center of mass system): 



ipis,{Q) = MH 



Mh Mh 



(20) 



where (?± = (0, 9), and Mh is the mass of the corresponding meson. The equations 



</'4(9)=V^,+(9)=0 



give the constraints on the components of the wave function: 
'^3{1) = : , Ml) = 



^i{q)MH{ujQ+uJq) 



Then we can rewrite the relativistic wave function of state as: 



7o^i(^) + v'2(?) - <ii.M~^)- ^'^'^ '^"'^ 



75 • (21) 



Prom this wave function we can obtain the wave functions corresponding to the positive and the negative 
projection, respectively: 



Mh 



^Q~^q \ f^Q-^q 



+ 7o 



Mh 



{rUqi^Q + m^Wg) 
-991(g) +(^2(9)- 



+^q 



75 



— L0„ \ I m„ — m, 



ra„ — ma / \ uj„ — uj, 



7o 



rua 



niqUQ + rriQUJq 



{22) 



{mqu; + m ujq) 



^mQ^-mq 

+^q 



75 



(23) 
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Table 2: Three sets (1-3) of input parameters. A is in the unit of GeV^, others are in the unit of GeV. 



Set 


a 




A 


^QCD 


rrib 




rus 


rrid 


ruu 


(1) 


0.060 


-0.60 


0.20 


0.26 


5.224 


1.7553 


0.487 


0.311 


0.305 


(2) 


0.055 


-0.40 


0.19 


0.24 


5.130 


1.660 


0.428 


0.285 


0.278 


(3) 


0.063 


-0.787 


0.21 


0.275 


5.3136 


1.845 


0.557 


0.352 


0.3465 



And there are two more equations from the reduced BS equation ()15() . which will give us coupled integral 
equations, and by solving them we obtain the numerical results for the mass and the wave function: 

dk 1 



{Mh -iVQ-U}2) 



(27r)3 2u)QUJq{EQmq + Egva^ 



X <^ {Eqimq + EqmQ){Vs - K,) 



Lpi{k){iOQLOq + m„mq - q )+ ip2{k)imqUJQ + m^ajg) 



-{Vs + K.) [Mk){mQ + mg){EQ + Eg) + v92(A;)K - w,)(^q - Eg) 

UJ„ -LOnI f dk 1 



M1)-V2{1)- 



m„ — TUn 



(27r)3 2ujQUjg{EQmg + Eg 



X \ {Egmq + EqmQ){Vs - K) 



ipi{k){ujQu;q + rriQiTiq - q )- (^2 ( ) (m^Wg + m^Wg) 



(24) 



{Vs + K) [ipi(k)im^ + mq){EQ + Eq) - (^2(1)^ - w<?)(£^Q - ^9)] 9 • 1} , 



where Eq - ^/ -r n.^ ^..^ ^q - y ...q . ..^ 



+ k'i . Finally the normalization condition is 



(2vr)^ 



:4^l{q)^2{(l)Mfj { 



mr^-rUa 2q {lOqUIq +U}qmqj 



+ 



UJq {ujQmg^ujqm^f 



2Mh 



(25) 



(26) 



4 Average Kinetic Energy of Heavy Quark inside Heavy Mesons in 
State 



The average kinetic energy of the heavy quark inside heavy meson in state, in the BS method, is 



proportional to the average spatial momentum squared: 

2 I dqq -> J a; -Wg m -mq 2q {uqUI + ujqiiiq) 

■ - 2(f^{q)(f2{i)MH ^ ' ' 



+ 



+ 



(27) 



(27r)3 ^ ' ^^^^ ' ^^'"^ \mQ-mq ' uj^-ujq ' {oj^mg + ujqm^)'^ 
In order to solve numerically the relativistic Salpeter equation, we use three different groups of input 
parameters (i.e. parameters for the potential and the masses of quarks), as shown in Table 2, from the 
best fit values [21]: 
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Table 3: Mass spectra and for heavy mesons in 0~ states with three sets (1-3) of input parameters. 
'Ex' means the results from experiments and 'ER' is the error of experimental values. 'Th' means 
the results from our theoretical estimate. 





Be 


Bs 


Bd 


Bu 




Ds 


Dd 


Du 


M GeV(Ex) 


6.4 


5.3696 


5.2794 


5.2790 


2.9797 


1.9685 


1.8693 


1.8645 


ER of Ex 


±0.4 


±0.0024 


±0.0005 


0.0005 


±0.0015 


±0.0006 


±0.0005 


±0.0005 


M GeV(Th)(l) 


6.296 


5.3654 


5.2804 


5.2778 


2.9791 


1.9688 


1.8687 


1.8655 


M GeV(Th)(2) 


6.304 


5.3670 


5.2804 


5.2762 


2.9795 


1.9691 


1.8699 


1.8650 


M GeV(Th)(3) 


6.292 


5.3656 


5.2806 


5.2788 


2.9799 


1.9690 


1.8673 


1.8650 


GeV^(l) 


0.955 


0.434 


0.348 


0.345 


0.615 


0.339 


0.280 


0.277 


/i^ GeV^(2) 


0.906 


0.429 


0.354 


0.350 


0.596 


0.331 


0.280 


0.277 


fi'i GeV^(3) 


0.958 


0.446 


0.350 


0.347 


0.636 


0.352 


0.286 


0.284 



a = e = 2.7183, a = 0.06 GeV, Vq = -0.60 GeV, A = 0.2 GeV^, Kqcd = 0.26 GeV and 
mfe = 5.224 GeV, rric = 1.7553 GeV, nis = 0.487 GeV, rud = 0.311 GeV, ruu = 0.305 GeV. 

With these three input parameter sets, we now solve the full Salpeter equation and obtain the masses 
and wave functions of the ground 0~ states. We list the calculated mass spectra of some 0^ states as 
well as the measured experimental values in Table 3. Then, by using the obtained wave function of 
heavy meson, we calculated /u^ from Eq. ()27() . as shown in Table 3. As can be seen from Table 3, if we 
change the values of the input parameters (sets 1-3) used in solving the Salpeter equation, we find that 
the obtained values of /i^ are almost unchanged (especially for Bd, Bu, Dd and Du mesons) when these 
parameters give a reasonably good fit of mass spectra. Therefore, we notice that our results for ^u^ are 
quite insensitive to the model parameters within the instantaneous BS method. We also note that the 
average kinetic energies of the heavy quark in different mesons differ significantly even when the heavy 
quark is the same, e.g. the value of ^-J^ of the heavy quark is significantly larger in Bg meson (~ 0.43 
GeV^) than in Bd (~ 0.35 GeV^) or Bu meson 0.34 GeV^). The difference of about 0.08 GeV^ is 
not a value which can be ignored compared with the value of /^^ itself. The bigger value of /i^ inside Bs 
meson than inside Bd or Bu means that b quark has a smaller residual momentum in Bd or Bu than in 
Bg. This implies that h quark is bounded more deeply in Bd or B^, than in Bg meson. In other words, 
the kinetic energy of the same h quark in heavy meson is more restrained by a light partner quark than 
by a heavy one, which is consistent with the running behavior of a^. Since our calculation of the average 
kinetic energy of the heavy quark has used the relativistic wave functions obtained from the full Salpeter 
equation, our results of the average kinetic energy are relativistic. Note that our results are quite 

9 



Table 4: The calculated uncertainties (in per cents) if we allow changes of all input parameters simulta- 
neously within 5% of the central values. 





Be 


Bs 


Bd 


Bu 


Vc 


Ds 


Dd 


Du 


AM/M 


±6.5 


±6.0 


±5.8 


±5.8 


±7.2 


±7.5 


±7.3 


±7.2 




±7.0 


±12.0 


±10.6 


±10.7 


±9.5 


±10.3 


±10.3 


±10.5 



different from the previously estimated ones of the potential model |13[ I14| I15j . This shows that the 
relativistic corrections are quite large, and cannot be ignored. 

In Table 4, we also show the calculated theoretical uncertainties for our results of the mass and average 
kinetic energy when we allow variations of all the input parameters simultaneously within 5% range of 
the central values. In comparison, our result for B^ d 

/x^ f» 0.31 - 0.38 GeV^ [our estimate] 

is a little larger than the recently experimentally derived CLEO values of 

/i^ = 0.25 ± 0.05 , [131 

and 

nl = 0.24 ±0.11 . 122 

In conclusion, we calculated the average kinetic energy of the heavy quark inside B or D meson by 
means of the instantaneous Bethe-Salpeter method. We solved the relativistic Salpeter equation and 
obtained the relativistic wave function and mass of 0~ state. Then we used the relativistic wave func- 
tion to calculate the average kinetic energy of the heavy quark inside the heavy 0~ state. We obtained 

(= -Ai) ^ 0.35 {B^,B^), 0.28 {D°,D^), 0.43 {Bs), 0.34 {Ds), 0.96 {Be) and 0.62 {rj^) GeV^. 
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Grant No. R02-2003-000-10050-0 from BRP of the KOSEF. The work of G.W. was supported in part by 
BK21 Program and in part by National Natural Science Foundation of China. 
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